On meson dominance in the 'second class' r — > rjirur decay 
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Abstract 



Motivated by recent estimates of the isospin-violating process r — > nirv T in QCD, mostly relying 
on the p and ao dominance of the relevant form factors near threshold, we present an assessment for 
the branching ratio that accounts for additional, potential, effects from the lowest radial excitations 
p' = p(1450) and a' Q = ao(1450), respectively, also lying in the decay phase space. 
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Renewed interest has recently been devoted to the semileptonic decay TjH- t]tiv t . This 
process belongs to the historical category of "second class current" decays [1], and is only 
viable in the Standard Model V-A semileptonic weak interactions through isotopic spin and 
G-parity violation, with a probability strongly suppressed by the smallness of the u — d 
current quark mass difference. Therefore, if measured with a significant branching ratio, 
this decay could signal some "new physics". Current theoretical estimates indicate for this 
transition a Standard Model branching ratio of the order of a unit times 10~ 5 2|-[5|, to be 
compared with the experimental upper limit 1.4 x 10 -4 [6|]. This should leave some room for 
revealing "non-standard physics" in high statistics studies of this decay, as recently discussed 
in Ref. 2| and, also, earlier in Ref. [7]. 

The theoretical estimates of the r —¥ r\'KV r form factors in the timelike region are based 
either directly on ad hoc p(770) and ao(980) polar forms or, like in Refs. [3|, [5(, on chiral 
perturbation theory [8| supplemented, above the t]tt threshold, by p and a poles. In this 
note, we present an attempt to phenomenologically parameterize the form factors in terms 
of both the chiral SU(3)xSU(3) symmetry breaking expansion and the p, ao resonance poles 
with, in addition, account of the lowest p, a® radial excitations p(1450) and ao(1450). These 
resonances are observed experimentally with masses well- within the rj7r invariant mass range 
allowed to the decay and, therefore, it should be interesting to assess the size of their effects 
relative to the ground state meson exchanges. Among other things, this would allow to make 
tests of the earlier theoretical predictions. 

With Vfj, = u'jfxd the weak vector current, the hadronic matrix element for r — y TyKv T can 
be decomposed into spin-1 and spin-0 exchange form factors, respectively, as 
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Here, q = p + k and t = q 2 is the invariant mass squared of the emitted rjTT pair. Finiteness 
of the matrix element at t = implies the condition /i(0) = /o(0). 

With t = (M v + M n ) 2 and X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + yz + zx), the partial width 

is: 
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For the value of form factors at t — 0, the prediction from leading and next-to-leading 



order in chiral symmetry breaking, explicit 
the process at the current quark level, is [3|, 



displaying the isospin breaking character of 
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In Eq. ([3]), e V7T is essentially the 77 — 71"° mixing angle, and fa denotes the average non-strange 

(2) 

current quark mass, fa = (rrid + ra„]/2. The next-order correction to t^-l is found to be 
reasonably small, e^il /e rj7 l — 0.25 [3j, |9J]. We include this correction at t = in Eq. (jHJ), and, 
consequently, will in the sequel adopt the value e vn = 1.34 x 10~ 2 , corresponding to the 
quark mass ratios ra u /raci ~ 0.55 and rag/nid ~ 18.9 [lOj. 1 

The information on the needed p' and a' coupling constants is scarce. For the radial 
excitation of the p, in practice only M p > ~ 1.46 GeV and V p i ~ 0.4 GeV are found in PDG 



ll| . We model the form factors for t ^ by expressions where the scale of the process is set 



by the isospin breaking parameter e V7T as an overall factor, times resonant 1~ and + poles 
evaluated in the SU(2) (or the SU(3)) limit. This kind of approximation will be used for the 
relevant coupling constants only, supposedly smooth, but not for the phase space in Eq. 

Accordingly, for the spin-1 form factor we assume, along the lines of Ref. 
tracted polar form with p and p' 
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where, denoting by M and V the resonance mass and total width, the threshold behaviour 
of p, p' — > rjii is attributed to T(t): 
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Here, L = 0, 1 for spin-0 and 1, respectively, and q is the momentum in the tjtt CM. frame. 
A similar parameterization was proposed in an approach bease on hadron duality and large 
number of colors iV c , in Ref. [13[. Also, in a sense, Eq. (T4|) resembles the modification of the 
p propagator introduced in Ref. [141 ] . 



Alternatively, on can consider a once-subtracted polar form, also satisfying /i(0) = e 
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1 Also, isospin breaking of electromagnetic origin will here be neglected. 



In the above equations, /'s and g's are the vector meson couplings to e + e~ and to rjir, 
respectively, for which we have made the assumption g p ( p ') v -K = e ring P ( P ')Tm- Numerically, this 
determination of the coupling constant g mn is qualitatively consistent with the estimates, 
obtained by different approaches, of Refs. 2| and [15l.ll6j]. The coefficient (3 P in front of the p' 
contribution is defined as (3 p = (M p /M p /) 2 x (f p ,g p , vlT f f p g p7Tn ). For a qualitative, theoretical, 
upper limit on the p' effects, a calculation within the constituent quark potential model 
gives fp'/fp — 1.10 from the ratio of the wave functions at the origin, see, e.g., Ref. |l7| . 
One can assess the value of g p ' nn (actually, an upper limit) by identifying T(p' — > tttt) to 
Y p i = 0.4 GeV, so that the upper limit on \i] p \ is derived: \(3 P \ < 0.18. One can notice that, 
for f3 p ~ —0.17, the analogues of Eqs. (jlj) and ([H]) reproduce the branching ratio of about 
25% for t — > ii7iu T 11], as well as the fit of the pion for factor from this decay 18j. Finally, 
from experimental data on p decays [ll], one obtains the ratio f p g p7T7T /M p ~ 1.2 needed in 
Eq. ([6]), with a small uncertainty. 

Apparently, of the two kinds of form factor parameterisation introduced above, Eq. (jl]) 
should be preferable, as involving a minimal number of input parameters and being less sen- 
sitive to the larger t-behaviour. Nevertheless, also the results from Eq. ([6]) will be presented, 
as potentially giving some useful information on the role of the form factor tails and (upper) 
limits on the effects due to radial excitations in the approach followed here. In Tab. [H we 
show the predictions for the spin-1 contributions to the decay branching fraction, for (3 P = 
(no p'), and for the maximal value \/3 p \ = 0.18 obtained above. 



TABLE I: Spin-1 branching ratio vs. (3 p . 
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The /3 = prediction from Eq. (jl]) is in the ballpark of the earlier estimates, and the 
addition of the p' could double the spin-1 branching fraction. Not unexpectedly, somewhat 
larger values, still of the same 10 -6 order of magnitude, are obtained from Eq. ([6]). One 
should notice from Eq. ((2]) that the contribution of /i(t) is suppressed, relative to fo(t), by 
the behaviour of phase space. 

We now turn to the spin-0 form factor /o(^) ; which we assume to be dominated by a 



and a' poles, similar to Eqs. fl4]) and ([6]). The nature of these scalar mesons is not well- 
established yet, and we will here assume them to be qq bound states. From Eq. ([I]), fo(t) is 
determined by the matrix element 

(« + (k) V (p)\ i0%|O) = (m d - m u )(7r + (k) V (p)\S^\0) = V2 (M 2 - M 2 ) f (t), (7) 

with S^ 1+l2 ^ = ud the quark bilinear scalar density. With the definitions 
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we modify the ao propagator to the unsubtracted polar form satisfying /o(0) = /i(0) 
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Here, the constant f ao is defined by F ao = e vn f ao explicitly showing the isospin violation of 
Eq. (jHJ), and T(t) is according to Eq. (jSJ). Moreover, g^rr* and g a ' ^ are the coupling constants 
for the decays of ao and a into 7/71", respectively, and f3 a is the ratio F a > ga'^/F^ g aovn . 

The coupling for a — > r\ix has been measured, with the result ga nmr — 2.80 GeV ac- 



companied by a total width T ao ~ 100 MeV [19], compatibly with [11]. As for the g a 



,r/7r 



coupling constant, we assume that the r/7r, t]'tt, KK and CU7T7T decay channels of the a satu- 
rate the total width T a > = 265 MeV. Using the ratios between partial widths r (7/71" )/r (7771"), 
r(KK)/T(j]7r) and r(u;7r7r)/r(?77T) reported in 111 ], one would find T(a — > t]ir) ~ 20 MeV, 



and the consequent estimate g a 'rm — 1-32 GeV. This value is in qualitative agreement with 



the estimate in Ref. 



20J derived from an SU(6) breaking scheme. 



To evaluate (3 a , we further need a theoretical assessment of the ratio F a > /F ao . One 
possibility could be to use, along the lines of Ref. [2j , current algebra equal-time commutators 
in an SU(6) framework to relate the + constant F ao to the analogous 1 + constant f ai , and for 
the latter the 2nd Weinberg sum rule f ai = f p . This would lead to — F ao ~ (md—m u )f p /M 2 , 



and a similar relation for F„' , so that: 
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By using the numerical value for the ratio f p >/f p suggested by the previous arguments made 
in connection to Eqs. (J3J) and ([6]), one would obtain the estimate (3 a ~ 0.23. 

An alternative is represented by the "linear dual model" ex plo ited in the framework of 
QCD finite-energy sum rules and local hadron duality in Ref. [2l| which, with S denoting 
+ 'scalar', indicates the rules: M^(n) = (n + 1)M| in perfect agreement for n = 1 with the 
measured M ao , and M a > Q , and fs( n ) = fs/( n + 1) hence F a ' o /F ao = 0.71. Numerically, this 
implies the estimate /3 a ~ 0.33, compatible within our approximations with the previous 
one. 

Finally, for the value of f ao , we take the QCD sum rule prediction F ao ~ 1.6 MeV 22]. 

These determinations enable us to numerically exploit Eqs. §§§ and (ITUj) in the assessment 
of the spin-0 part of the branching ratio of interest here. The results, for (3 a = and 
\/3 a \ = 0.33 (the largest value estimated above), are displayed in Tab. [TTJ 



TABLE II: Spin-0 branching ratio vs. /3 a . 
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The size of the spin-0 branching ratio can be one order of magnitude larger (0(1O~ 5 )) 
than the spin-1 part. For f3 a = (no a' ), the estimated value is qualitatively consistent with 
earlier estimates. Including the a (1450) contribution, and depending on the value of /3, the 
rate might be increased by a factor of two. 

In conclusion, we confirm the orders of magnitude of 10~ 6 and 10 -5 for the spin-1 and 
spin-0 branching ratios of previous Standard Model estimates, these predictions therefore 
seem rather robust. However, the effects from the radial excitations p(1450) and ao(1450) 
may be non-negligible and amount to a factor of two or so, compared to the simpler p and 
ao polar forms. The upper limit from the above tables is BR(r — > T]7rv T ) ~ 2.6 x 10~ 5 . 
On the one side, this value could be in the reach of future high statistics experiments on 
r-decays. On the other, considering the current experimental upper limit of order (9(10~ 4 ), 
there could still remain some room for an unexpected discovery. 
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